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Abstract 

We study tunneling of weakly-interacting spinless electrons at zero temperature through a single <5 barrier in one-dimensional 
wires and rings of finite lengths. Our numerical calculations are based on the self-consistent Hartree-Fock approximation, 
nevertheless, our results exhibit features known from correlated many-body models. In particular, the transmission in a 
wire of length L at the Fermi level is proportional to L~ 2a with the universal power a (depending on the electron-electron 
interaction only, not on the strength of the 8 barrier). Similarly, the persistent current in a ring of the circumference L obeys 
the rule / oc L~ 1— a known from the Luttinger liquid and Hubbard models. We show that the transmission at the Fermi 
level in the wire is related to the persistent current in the ring at the magnetic flux h/4e. 
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It is known that a one- dimensional (ID) wire which 
is free of impurities and biased by macroscopic contacts 
yields the quantized conductance. This effect can be 
derived within a simple model of non-interacting elec- 
trons if we assume a negligible electron back-scattering. 
Placing a single impurity (scatterer) into the wire, the 
conductance is not quantized any more. It is because 
of the back-scattering of the electrons from the scat- 
terer. The wire conductance (given by the Landauer 
formula) is proportional to the electron transmission 
at the Fermi level [1]. 

When the electron-electron (e-e) interaction is con- 
sidered, the conductance of an infinitely long wire with 
a single scatterer decreases as T~ 2a for (T — ► 0). This 
behavior is known from the Luttinger liquid model [2] 
where the power a depends on the e-e interaction only. 



Assuming a repulsive interaction (a > 0), the scatterer 
becomes impenetrable at zero temperature regardless 
the strength of the scatterer. 

Matveev et al. [3] studied the Landauer conductance 
of the interacting ID electrons through a S barrier 
in a wire with contacts. They analyzed the effect of 
the Hartree-Fock potential on the tunneling transmis- 
sion assuming a weak e-e interaction. They derived the 
transmission using the renormalization group (RG) ap- 
proach and confirmed the universal power law T~ 2a . It 
is believed that this approach goes beyond the Hartree- 
Fock approximation. 

In this paper we consider the non-Luttinger liquid 
model of the same type as analyzed by Matveev et 
al. [3]. However, we do not use the RG theory. We apply 
the self-consistent Hartree-Fock solution by means of 
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numerical calculations instead. We calculate the trans- 
mission probability at zero temperature which is re- 
lated to the Landauer conductance. A good agreement 
with the theory of Matveev et al. is found. In particu- 
lar, we simulate dependence of the conductance on the 
wire length (L) for various 5 barriers and we reproduce 
the universal power law oc L~ 2a that becomes asymp- 
totic for large wire lengths and/or strong S barriers. 

We also consider a mesoscopic ring threaded by the 
magnetic flux. As a consequence, the persistent cur- 
rent (7) arises [4]. We study the persistent current of 
interacting spinless electrons with a single S barrier at 
zero temperature. We show how the transmission can 
be extracted from the persistent current in the ID ring 
where / oc L~ a ~ 1 and compare it with the transmis- 
sion obtained from the ID wire. 

First, consider ID wire of length L with interacting 
spinless electrons. Both wire ends are connected to con- 
tacts and the single S barrier is localized in the center 
of the wire. The single-electron wave functions ip k (x), 
where k is the electron wave vector, are described by 
the Hartree-Fock equation Htpk(x) = ekipk(x). The 
Hamiltonian has the form 

H = "£ ^+-yS(x-L/2) + U H (x)+U F (k,x), (1) 

where ^5(x — L/2) represents the localized scatterer. 
We apply the boundary conditions 

»/. ( ~ n\ ikx | —ikx 

ip k (x = 0)=e +r k e 
ipk{x = L) = t k e kx , 
ip-k{x = 0) = t' k e~ lkx , 

^- k (x = L)=e- ikx +r' k e ikx (2) 

with r k and t k being the reflection and transmission 
amplitudes, respectively and k > 0. Knowing the trans- 
mission, we can easily evaluate the Landauer conduc- 
tance {e 2 /h) \tk F | ■ 

The Hartree potential induced by the S barrier reads 

L 

Uh{x)= J dx'V(x-x') 
o 

k p 

x / ^[\^{x')\ 2 -WUx')\ 2 ] 0) 

— kp 

and the Fock non-local exchange potential is written as 



L 

Up(k,x) = — r~rr dx'V(x — x') 
Vk{x) J 
o 

k p 

f rlk' 

x J ^Mx'M>(x')1>*(x) (4) 
-k F 

with V(x — x') being the e-e interaction. Essentially 
the same one-dimensional model was considered by 
Matveev et al. [3] . 

The wire is connected to large contacts via adia- 
batically tapered non-reflecting connectors [1]. In the 
case of interacting electrons without presence of the 
scatterer (7 = 0), we assume that there is no back- 
scattering at the wire ends due to the adiabatically 
tapered connectors. Then, the solution of Eq. (1) is 
described by the free wave, ipk(x) = e %kx , having the 
eigenenergy 

e k = h 2 k 2 /2m + U { F(k) (5) 

with Up(k) = UF[ipk{x) = ipk(x)] being the nonzero 
Fock shift. Note that this solution is valid if we implic- 
itly assume that the Fock interaction is present also 
in the contacts. If the energy in Eq. (5) holds inside 
the wire and we turn off the Fock shift to zero out- 
side the wire, we obtain at each wire end the potential 
drop U F (k). This would cause back-scattering at both 
wire ends and the solutions e lkx and e ~ lkx would be no 
longer valid (in contrast to the ballistic conductance of 
clean wires [1]). 

If the barrier yS(x — L/2) is positioned in the 
wire, this 5 barrier induces Friedel oscillations of the 
Hartree-Fock potential. The Friedel oscillations pene- 
trate through the wire ends into the contacts, where 
they decay fast due to the higher dimensionality of the 
contacts and decoherence. To mimic this decay within 
our model, we sharply turn off the oscillations to zero 
at both wire ends keeping Uh = and Uf — U F (k) 
outside the wire. Such a constant potential emulates 
the non-reflecting connectors and justifies the above 
boundary conditions. 

Our numerical results are carried out for the GaAs 
wire with the corresponding effective electron mass 
m = 0.067 mo, the electron density n — 5 x 10 7 m _1 , 
and the short range e-e interaction 

V(x-x') = V e-\ x - x '\ /d . (6) 

We use the short range e-e interaction (6) because of 
comparison with the RG theory [3] where the e-e in- 
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teraction is assumed to be finite. Physical meaning of 
the finite range is the screening. 

The asymptotic formula for the transmission prob- 
ability at the Fermi level derived for weak e-e interac- 
tions [3] reads 



\tk f 



\tk 



{d/Lf 



\~tkp 



\r kF r + \t kF \ (d/L)** \r kf 

(7) 

where d is the range of the e-e interaction V(x-x') and 
ik and f k describes the transmission and the reflection 
amplitudes of the bare S barrier [3], respectively. The 
right hand side of Eq. (7) remains valid for small tk F 
and/or large L. For weak e-e interaction (a <JC 1), the 
power a reads [3] 

a= V(0)-V(Zk F ) 
2nhvF 

where V(q) is the Fourier transform of the e-e interac- 
tion V(x— x'). We evaluate a for our e-e interaction (6), 
for which V{q) = 2V d/{l + q 2 d 2 ). 

The bare amplitudes are t k = k/(k + i£) and f k = 
—iC,/(k + where £ = jra/h 2 . As k F and m remain 
fixed, we parametrize the bare 8 barrier by its trans- 
mission coefficient \tk F j in the following. 

Second, consider a ID wire for which the Hartree- 
Fock equation reads 



2m 



.8 + 2ne 
dx Lh 



+ U H (x) + U F (j,x) 



<t>\ + jS(x-L/2) 



i>j(. x ) = £ji>j{ x ) (9) 



satisfying the boundary condition ij)j(x + L) = ipj(x). 
The persistent current is given as I = — gr-E(0), where 



£j - g i^i \Uh{x) + U F (j,x)\ tpj) 



(10) 



In case of the non-interacting electrons, the persis- 
tent current can be evaluated [5] as follows 



(ev F /2L)\t £F | sin(27re0/7i), 



(11) 



with the transmission amplitude of the scatterer at the 
Fermi energy \t SF \ -C 1 and the Fermi velocity v F . 
Assume L — > oo, the repulsive e-e interaction a > 
0, and <f> = h/Ae. Then, replacing \t SF \ by \tk F \ and 
inserting Eq. (7) into Eq. (11), we obtain [5] 
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Fig. 1. Conductance in unit of the transmission probability 

I - I 2 

at the Fermi level vs length for various <5 barriers \tk F \ 
in a log-log scale. The dashed curves show the data corre- 
sponding to the formula (7). The rilled circles connected by 
the full lines arc our self-consistent Hartrcc-Fock data for 
the wire; the curves approach the same asymptotic power 
law as the dashed lines. The open circles are extracted from 
the persistent current using Eq. (12). 



We solve Eqs. (1) and (9) by means of the self- 
consistent iterative procedure and follow approxima- 
tion in Ref. [6] in order to decrease computational time 
and memory. One can further simplify Eq. (4) to 



U F {x) 



■jdx'V(x-x) J ^Re[i>l,{x')i}y{x)] 

(13) 



-k F 



noticing that /_£ dk'ipy (x') ipk' ( x ) — 2tv8(x — x'). 
The Fock potential (13) becomes local and independent 
on k. This allows us to simulate longer wires than in the 
case of the non-local Hartree-Fock model, cf. Eq. (4). 
We also present the results without this simplification 
but for a substantially shorter wire lengths. 

Figure 1 shows the transmission probability |ife F | 2 
for the wire and the ring versus the length L for var- 
ious 8 barriers. The result of the RG (7) is presented 
by the dashed lines. For strong <5 barriers the dashed 
lines follow the asymptotic power law \tk F | 2 oc L~ 2a as 
manifested by the linear decay with slope —2a in the 
log-log scale. Our Hartree-Fock curves (filled and open 
circles) show slightly higher transmission following the 

I ~ ] 2 

same slope —2a. Note that for | tk F | small enough and 
substantially longer wires and rings, all the Hartree- 
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Fig. 2. A log-log plot of the persistent current vs ring cir- 
cumference (left) and the conductance (transmission) vs 
the wire length for different a. 



Fock curves decay with the same slope as the dashes 
curves, independently on the strength of the 8 barrier. 

In Fig. 2 we present the squared persistent current 
(which has the meaning of transmission, cf. Eq. (12)) 
and the conductance as functions of the length L for 
both the ring (left) and the wire (right). The data 
shown were calculated for the bear 8 barrier \t kp | = 
0.003 and for various e-e interactions. As seen, both 
data for the wire and the ring follow the same behavior 
having the slope —2a. We have chosen the e-e interac- 
tions a = 0.0277, 0.0561, and 0.0855 that correspond 
to Vo = 11 meV, 22.3 meV, and 34 meV, respectively, 
and for the common range of e-e interaction d — 3 nm. 

The presented results were calculated within the lo- 
cal Fock term (Eq. (13)). When the non-local Fock term 
is considered, the resulting effective potential gains a 
nonzero imaginary part which changes the transmis- 
sion properties of the system studied. In Fig. 3 we de- 
picted the conductance at short lengths for the wire 
(triangles) and the ring (circles) within the local ap- 
proximation [6], the wire for the non-local Fock term 
(squares) and the non-local Fock term for which the 
imaginary part is omitted (stars). 

If the Fock term is completely omitted, we obtain 
different results for the conductance (or the persistent 
current). It is because the perfect reflection is resulted 
instead (in the case L — ► oo). However, the power-law 
behavior is expected to persist if the density functional 
theory is applied. 

Note that we have not yet reached the asymptotic 
regime because the conductance in Fig. 3 is not linear in 
the log- log scale. Longer wire lengths are required when 
the non-local Fock term is involved in our calculations. 
Unfortunately, the memory limitations (4 GB) do not 
enable us to obtain reliable data for longer systems. 

In conclusion, we have calculated the Landauer con- 
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Fig. 3. Conductance vs. length of the system for the wires 
and the rings within local and non-local Fock approxima- 
tions for a = 0.0855 and \t kF = 0.005. The dashed line 
shows the analytical formula (7) of Matveev et al. 

ductance and the persistent current of the weakly- 
interacting spinless electrons in the ID wires and rings 
with a single 8 barrier. We have used the self-consistent 
Hartree-Fock approximation at zero temperature. We 
have found the universal power law \t kF \ 2 <x (d/L) 2a 
known from the Luttinger-liquid model [2] and the RG 
models [3,7]. We conclude that the universal power law 
is not exclusively the consequence of correlations; it 
can be obtained within the Hartree-Fock approxima- 
tion. To prove this, longer lengths are needed. We have 
also found that essentially the same wire conductance 
can be extracted from the persistent current in the ID 
rings. 
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